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DEVELOPMENT OF A TRANSFER FUNCTION METHOD 


FOR DYNAMIC STABILITY MEASUREMENT 

Wayne Johnson 

Ames Research Center 
and 

Ames Directorate, USAAMRDL 


SUMMARY 


A flutter testing method based on transfer function measurements is 
developed. The error statistics of several dynamic stability measurement 
methods are reviewed. It is shown that the transfer function measurement con- 
trols the error level by averaging the data and correlating the input and out- 
put. This method also gives the experimenter a direct estimate of the error 
in the response measurement. An algorithm is developed for obtaining the 
natural frequency and damping ratio of low-damped modes of the system, using 
integrals of the transfer function in the vicinity of a resonant peak. Guide- 
lines are given for selecting the parameters in the transfer function measure- 
ment. Finally, the dynamic stability measurement technique is applied to data 
from a wind-tunnel test of a proprotor and wing model. 


INTRODUCTION 


A major task in developing new aircraft is the demonstration of satisfac- 
tory aeroelastic characteristics. Beyond the determination of the flutter 
boundary, or establishing that the aircraft is flutter-free throughout its 
operating envelope, it is also desirable to obtain detailed dynamic informa- 
tion. This information is used to determine specific characteristics (such as 
the aircraft gust response) and to verify mathematical models of the aircraft. 
Thus, the flutter testing of an airplane or helicopter requires an accurate, 
efficient, and reliable method to measure the aeroelastic response, and then a 
method to obtain from the response the parameters defining the dynamic 
characteristics. The most important parameters are generally those defining 
the flutter stability level, that is, the frequency and damping of the low- 
damped modes of the system. The state of the art of dynamic stability testing 
in the aircraft industry has been described in a number of surveys (refs. 1 to 
5) . The problem usually involved is the flutter testing of the airplane wing 
or tail in flight. Some form of frequency-response procedure is common. Many 
variations exist, but a typical procedure involves fast swept-sine excitation, 
digital analysis of the response (analog analysis is also still common) , with 
the frequency and damping determined from the circle on the phase plane (see 
ref. 6). Another approach often used is to determine the dynamic stability 
from the decaying transient response. A number of measurement and analysis 
techniques has been developed for structural dynamics testing, but in such 
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testing the principal data required are the frequencies and mode shapes, 
whereas in flutter testing it is the damping that is of primary concern. 

Flutter testing methods such as fast sine sweep or transient decay can 
work well with clean data. With low process and measurement noise, the excita- 
tion level can be sufficiently high that the effects of noise on the measure- 
ments can simply be neglected. As the noise increases, however, the accuracy 
of the flutter parameter determination by such methods degrades rapidly, espe- 
cially the damping measurement. Dynamic stability measurement is particularly 
difficult for rotorcraft, where many degrees of freedom are involved, and the 
levels of process noise (turbulence- induced , especially in a wind tunnel) and 
measurement noise (rotor- or engine-produced vibration) can be high. Thus, in 
general, a dynamic stability measurement technique is required that can be 
applied to data with a significant noise level. 

Accurate and reliable flutter testing requires minimal error in measuring 
system response and in analyzing the data. Regarding control of the error 
level in the response estimation, the test techniques may be classified in 
three categories. First are the techniques that use a high signal-to-noise 
ratio to control the error. Examples are the transient decay method (perhaps 
using the moving block analysis in the data reduction) , and a fast sine-sweep 
measurement of the transfer function. Second are the techniques that average 
the estimation of the response due to existing disturbances of the system 
(such as aerodynamic turbulence). Examples are spectral analysis, correlation, 
and random decrement signatures. Third are the techniques that, in addition 
to averaging, use correlation between the response and a measured external 
input to reduce the error. Examples are single-input, single-output transfer 
function methods, and more sophisticated parameter identification techniques 
using several output and input measurements. The error characteristics pro- 
gressively improve with each of the three categories. After the response meas- 
urement is obtained, in either the time or frequency domain, the parameters 
describing the dynamic characteristics of the system must be calculated, par- 
ticularly the damping ratio, which gives the quantitative level of stability. 

A technique is required that is accurate when applied to noisy data and that 
can be implemented for online stability measurements. Most techniques 
(including the one developed here) obtain the frequency and damping of the 
low-damped modes by an algorithm based on the response of a single-degree-of- 
freedom, second-order system. 

In this report, the method for measuring dynamic stability uses a single- 
input, single-output transfer function, and pays particular attention to the 
error characteristics with noisy data. The error statistics of several 
response measurement techniques are summarized to establish the relative value 
of the transfer function. Then, an algorithm is derived for obtaining the 
damping ratio and other parameters. Next, the procedure for measuring the 
transfer function is discussed. Finally, an example of the application of 
this dynamic stability measurement technique is given. The data used are from 
a wind-tunnel test of a model proprotor and cantilever wing. We begin with a 
discussion of the mathematical model used for the aeroelastic system. 
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LINEAR SYSTEM DYNAMICS 


A linear system is considered in the examination of the dynamic stability 
measurement techniques, with excitation by various control inputs and external 
disturbances. The response of the system is measured, perhaps with signifi- 
cant measurement noise such as rotor- or engine- induced vibration. The system 
can be excited by existing unknown disturbances (such as aerodynamic turbu- 
lence) , or by a measurable external input applied to determine the dynamic 
stability. The system motion is therefore described by linear, time- invariant 
differential equations, of the form 

'k = Ax + Bu 

y = Cx + v 

where x is the state vector (degrees of freedom), and ^ is the measurement. 
The vector ^ is random measurement noise; is an input exciting the system, 
either an existing random disturbance or an external input (random or deter- 
ministic). The matrices A, B, and C are constant since the system is time 
invariant. The solution for the response to excitation by ii, with initial 
conditions at time to* is 

*(t) - e A ( t - t o) J( t 0 ) + f' e A(t - T> BS dT 

t 0 

jKt) = C$(t) + V 

(see ref. 7). The stability of the system is determined by the eigenvalues 
of A. The eigenvalues usually occur in complex conjugate pairs, of the form 

A = -£27TO) n ± i27T(jo n /l - C 2 

where u) n is the natural frequency (in Hertz) and £ the damping ratio of 
the mode. The mode is stable if £ > 0, and is an exponentially decaying 
oscillation for 0 < £ < 1. In flutter testing, the low-damped modes of the 
system are of primary concern (roughly 0 < £ < 0.1). If A is the diagonal 
matrix of the eigenvalues of A, and M the matrix whose columns are the 
corresponding eigenvectors, then A = MAM -1 . For further information on the 
dynamics of linear systems, the reader is directed to references 7 to 9. 

Regardless of the procedure for measuring stability, the system must be 
controllable and observable. Thus, the input must sufficiently excite the 
modes of interest and the modes must be observable in the response of the 
variables measured. 
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ERROR ANALYSIS OF THE RESPONSE MEASUREMENT 


Error Analysis of Transient Decay Method 

The transient motion of a dynamic system is composed of exponentially 
decaying oscillations of each mode. A particular measurement is often domi- 
nated by one low-damped mode. From the oscillation period and the decay rate 
of the trace, the frequency and damping of the mode may be estimated. The 
typical procedure establishes a large sinusoidal motion by external excitation 
at the natural frequency of the mode of interest. The excitation is stopped 
and the subsequent transient motion is analyzed to determine the damping. 

Consider the linear, time- invariant system described in the previous sec- 
tion, with state variable i and measurement y; u is process noise, such as 
aerodynamic turbulence; and v is measurement noise. It is assumed that $ 
and v are random disturbances with zero mean. The motion following time tQ , 
where the initial conditions x(to) are established by the external excitation, 
is 

f(t) = Ce A(t_to) x(t 0 ) + C / 1 e A( ' t ~ T - > Bi dx + v 

The first term in y is the transient we wish to observe. The second term is 
the noise in the response caused by the disturbances occurring after tQ. The 
third term is the measurement noise. The error introduced by these disturb- 
ances is the primary difficulty with this technique. The following analysis 
derives the expected value and variance of the response, designated Ey and 
Vy = E (y - Ey) 2 . The terms "expected value" and "variance" are used rather 
than "mean" and "standard deviation," because the latter are often associated 
with a time average of the response. Here, Ey and Vy are themselves time- 
varying quantities. The normalized variance e = VVy/Ey is the basic error 
parameter considered. For a discussion of the probabilistic aspects of such 
an error analysis, see reference 10. 

The expected value of the measurement is the transient response 

Ey = Ce^^ to ^x(to). The external excitation must establish a large, nonzero 
initial condition for the modes of interest. The variance of the observation 
is : 


Vy = E (y - Ey) 2 = E(^ 2 ) - (Ey) 2 

= c/ t e A( ' t-T ' ) BQB T e AT( - t-T LT C T + R 
tn 


where it is assumed that u is white noise with correlation 

Eu(t 1 )u t (t 2 ) = Q<Ht! - t 2 ) ( so also u for t > to is uncorrelated with x 

at tg), and the measurement noise covariance is R = E [v(x) v^ ( t) ] . The steady- 
state variance X = Vx of the response of the system to the disturbance u 
alone is related to Q by 

-(AX + XA t ) = BQB T 
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(ref. 7). If this relation is substituted for BQB T in Vy, then the output 
variance is proportional to the square of the standard deviation of the steady- 
state response cr x 1 2 = X. In particular, if the response is dominated by a 
single mode, then the normalized variance is approximately 

e 2 = JL 4 /^\ 2 r e - 2 Mt-t 0 ) _ 1 1 + r_ e - 2 A(t-t 0 ) 

(Ey ) 2 \ x 0 / L J x 0 2 c 2 


(see ref. 11). The first term is the process noise; the second term is the 
measurement noise. As t - tg increases and the transient decays, the noise 
becomes more important. To determine the damping ratio, the transient must be 
observed while it decays to a fraction f of the initial value (f = 0.3 to 
0.5, typically). The maximum error to be dealt with is thus 


2 . /£x\ 2 

e max - 


(f -2 - 1 ) + 


2 2 
*0 c 


f-2 


which is basically the ratio of the background process and measurement noise, 
and the smallest value of the transient. The maximum rms error is then at 
least two to three times the ratio of the rms response due to u and the 
initial response amplitude, plus the measurement noise. 


The frequency and damping from the transient decay are usually based on 
the response of a single mode 

x = xq e ^ 27 Ta 3 ti t cos ^ 2 iT(jL) n /l - £ 2 t + 4 > o ) 


The frequency is given by the period of the oscillation co = 1/T (Hz) . To 
determine the damping ratio, consider two peak amplitude measurements Xj an d 
X 2 , which are n oscillations apart. Then, since x^ = exp (- 27m c/ - C z ) > 
it follows 


C 


1 + 


27m 


in X1/X2 



£n x^/x2 
27rn 


Often, this task is accomplished by hand, directly from an oscillograph trace 
of the response. In general, the damping ratio may be determined from the 
envelope of the oscillation by £ = -(d £n x/dt)/2Traj. The difficulty in 
reducing the data lies in determining the envelope of the decaying oscillation. 


Another means of analyzing the decaying transient response is the moving 
block method, which extracts the envelope of the decay. A Fourier transform 
is applied to a block of data y(t) from t = x to t=x+T; the magnitude 
of the spectrum line at the natural frequency oo n , plotted as a function of x, 
gives the envelope and hence the damping ratio. Specifically, the operator 


1 P x-KT , x -i2irv(t-x) 

-a Jr <■ • - )e 


dt 
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is applied to the response, yielding the function y(x,v) = Ly(t) . The 
magnitude of the expected value of y at the frequency v = a) n is propro- 
tional to the exponentially decaying envelope 

I Ey v=0J n I ^ constant * e ^ 27ra3 n( T t o) 

(see, e.g., refs. 5 or 11). The moving block analysis can also use analog 
components (a tracking filter locked to the frequency of the desired mode, an 
rms integrator, and a log converter). The slope of £n|y| as a function of 
x gives the damping ratio of the dominant mode in the response. The princi- 
pal advantage of the moving block analysis is that it extracts the frequency 
and exponential envelope from the decaying oscillation, which greatly facili- 
tates the mechanization of the data processing task. Working in the frequency 
domain also reduces the effect of measurement noise, since only noise around 
the natural frequency oj n is important. However, the moving block analysis 
involves a linear operator, which can have no fundamental influence on the 
normalized error (ref. 11). The response of the system to u is mainly a 
superposition of oscillations at the frequency m n , which is transmitted 
through the tracking filter along with the transient oscillation, so the 
moving block analysis has little effect on the process noise. 

Another approach is to curve-fit the measured response to an exponen- 
tially decaying oscillation in the time domain. An iterative least-squared- 
error technique is usually used, and a multimode curve fit may be implemented 
if necessary. Convergence problems are very likely with a high noise level, 
however . 

A good estimate of the damping can be obtained by the transient decay 
method with clean data. The moving block or tracking filter analysis helps by 
mechanizing the extraction of the damping ratio from the data, and by reducing 
somewhat the measurement noise due to the filtering. The accuracy of the 
damping estimate degrades rapidly as the noise increases, however. With sig- 
nificant noise, it may be impossible to obtain any meaningful results. Often, 
the experimenter can do little to control the relative noise in the response 
measurement either. The error analysis above shows the only options are to 
increase the amplitude of the initial conditions or to reduce the process and 
measurement noise. Sometimes the noise can be influenced by the experimenter 
(e.g., by improving the rotor track or engine isolation to reduce the vibra- 
tion or by conducting a flight test only under calm atmospheric conditions to 
minimize the process noise). In a wind tunnel, however, little can be done to 
reduce the turbulence, which is likely to be high at maximum speed. A sub- 
stantial increase in the initial excitation is usually not practical either, 
because of structural or amplitude limitations in both the response and input. 
Clearly, to handle noisy data adequately, a different test technique is 
needed . 
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Error Analysis of Output Autospectrum Method 

To improve the accuracy of the system response measurement, the data must 
be averaged. One means by which such averaging can be introduced is to meas- 
ure the response autospectrum. Consider a system excited by a random disturb- 
ance, such as aerodynamic turbulence, which is not measured. Let S^) (cj) be 

the autospectrum of the kth measured time series of the response y. Then, an 
estimator of the true spectrum of the response is given by the ensemble 
average over K calculations of the spectrum: 

K 

M«> 

k=l 

(ref. 12). The spectrum S can be obtained by digital or analog processing. 
With a digital processor and hardware FFT units, spectrum analysis is more 
convenient than autocorrelation analysis; also, working in the frequency 
domain is better suited for data reduction. 

Generally, the excitation of the system consists of multiple inputs (e.g., 
several gust components) , so the output spectrum is related to the input 
spectra by 

S y «]£ |Hii 2 S v . 
i 

where H-j_ is the transfer function between the response y and the ith input, 
and S v ^ is the (unknown) input autospectrum. For the system considered, Hj_ 

is a rational function, with the denominator D(co) - the characteristic equa- 
tion - the same for all inputs. Thus, with Hi = Ni/D, 



i 


The poles of the characteristic equation, which give the modal frequency and 
damping, are of primary concern. Sy is seen to be a direct measure of D(oo). 
The spectrum Sy depends also on tne characteristics of the input spectra. 

To obtain the properties of the system alone, the input spectra must be reason- 
ably flat. This criterion is not too severe, since it need be applied only in 
the vicinity of the resonant peaks, which are very sharp for low-damped modes. 

A 

An analysis of the statistics of the estimator Sy of the autospectrum 
is available in the literature (see refs. 12 and 13). The expected value and 
variance of S are 

ES(to) s S (u>) + S" (to) 
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VS (to) = 

where Aw is the frequency resolution (in Hz) of the spectrum, and K the 
number of averages. 1 These results, which are essentially independent of the 
statistical properties of the signal, are applicable to general random dis- 
turbances. The normalized bias error is 

= Aa) 2 S" 
e b 24S 

The most critical case is at a resonant peak, where the greatest accuracy is 
required and the curvature S M is largest. Assuming that, in the vicinity of 
the peak, the spectrum may be approximated by a single-degree-of-f reedom sys- 
tem, one obtains S/S" = (l/2)£ 2 w n 2 . So, for a given bias error, the follow- 
ing frequency resolution is required: 

Aw = /l2e^ £w n 

Using a frequency resolution of Aco - (1/2) £w n gives = 0.02, which is a 
negligible bias error compared to the usual variance levels (twice that reso- 
lution is often satisfactory). The half-power bandwidth of the peak is 
Ao )2 = 2£w n , so this resolution corresponds to covering the bandwidth of the 
peak with four or five spectral lines. The sampling time for a single spec- 
trum measurement is related to the frequency resolution (in Hz) by T q = 1/Aw 
(ref. 12). Hence, the spectrum frequency resolution is an important test 
parameter, affecting both the bias error and test time. 

The bias error is easily made negligible by a proper choice of the 
frequency resolution. The remaining error is the normalized variance: 
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This error is inversely proportional to the square root of the number of 
averages, a standard result for sample means. The estimate of the spectrum 
may thus be made as accurate as desired by increasing the number of averages. 
The total record length required for the spectral analysis is the product of 
the number of averages and the sample time for a single spectrum measurement, 
which gives 

T _ JL = 1 

Aw Awe 2 


2 If N samples are collected at rate r, then the time required to 
obtain a single record is Tq = N/r. The maximum frequency in the spectrum is 
(Omax = r/2 (the Nyquist frequency) and the number of spectral lines is N/2, 
so the frequency resolution of the spectrum is Aw = r/N = 1 /Tq. Then, the 
total sampling time with K averages is T = KTq = K/Aw, giving K = TAw, or 
K = TB e in the notation of reference 12. 
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With the above result for the frequency resolution, T = 2/ (£a) n e: 2 ) ; the total 
time over which data are collected must be increased as the damping, frequency, 
or error decreases. This sample time expression is a fundamental result for 
the amount of data required to estimate the statistics of a random signal. It 
is thus applicable to all no-input dynamic stability measurement techniques, 
such as correlation (which is simply a transform of spectral analysis to the 
time domain) or random decrement analysis (discussed below) . 

To determine the system properties from the response spectrum, the input 
spectrum is generally assumed to be flat in the vicinity of each resonance, so 
Sy = |h| 2 S V = |h| 2 . For the low-damped modes that are of most interest, this 
assumption is required only over the narrow frequency range Au)/to n = 2£. 
Occasionally, the response spectrum can be corrected also for some known varia- 
tions in the input spectrum. Data reduction then proceeds by treating S y as 
the magnitude of the transfer function. For example, if Aa)2 is the half- 
power bandwidth of the resonant peak (where the autospectrum Sy = [ H [ 2 has 
fallen to one-half the peak value) , then the damping ratio can be estimated 
from 5 = A(jL)2/2a) n . Alternatively, the damping ratio can be calculated by 
integrating the spectrum through the resonance 

i A” s y “ 2 dt0 

TTW n 3 (SyOpeak 

(This expression can be used to estimate the damping of multimode systems as 
well by limiting the integration to a narrow range about each resonant peak.) 

Spectral analysis, by introducing averaging of the data, gives the 
experimenter control over the error level in the response measurement. How- 
ever, techniques that measure only the output have two major difficulties. 
First, as seen above, the output spectrum depends not just on the system but 
also on the input characteristics, so there is always some uncertainty whether 
the system parameters are being correctly estimated. Second, the total record 
time to measure the spectrum accurately is often impractically long, particu- 
larly for full-scale aircraft (which have low fundamental structural fre- 
quencies) . Consider, for example, a case with to n = 5 Hz, £ = 0.02, and 
£ ~ 0.2; the above expression gives T = 500 sec. 

An autocorrelation analysis is entirely equivalent to an autospectrum 
analysis, the former working in the time domain and the latter in the fre- 
quency domain. The statistics of an autocorrelation measurement are thus sim- 
ilar to the spectrum measurement statistics summarized above (see ref. 13). 

The method of random decrement signatures (randomdec) is another procedure 
for analyzing the response of a linear system to unmeasured random distur- 
bances, designed to estimate the transient response of the system. The error 
statistics of the random decrement method are basically the same as for 
spectral analysis (see appendix). Spectral analysis, correlation, and random 
decrement techniques are all members of the general class of methods for 
analyzing the properties of a system using the response to an unknown random 
disturbance, with averaging to control the error level. A basic character- 
istic of such methods is that a large amount of data is required to extract 
information accurately from the response signal; thus, for low-frequency 
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systems, the time to collect the data is very long. The examination of the 
autospectrum measurement method has shown the need for a procedure that is 
more efficient in reducing the error and that measures the characteristics of 
the aeroelastic system alone. 


Error Analysis of Transfer Function Method 

The accuracy of the measurement of the system response is further 
improved by measuring the input as well as the output. Beyond the influence 
of averaging the data, the error is also reduced by the correlation between 
the input and output; furthermore, the properties of the system alone are 
obtained, in the form of the transfer function between the output and the 
input. Usually, a measurable external input is applied specifically to excite 
the system for the flutter testing. The cross spectrum between the input and 
output is found, in addition to the autospectra. The ratio of the cross 
spectrum and the input autospectrum then gives the transfer function H(w) , 
from which the system parameters may be found. 

Consider again a linear system, but now with measurements of both the 
input and the output: 

t = a£ + b 3 + ft 

y = Cx + vq 

z" = U + V-J_ 

Here, u is the random input to the system (which is measured), w is process 
noise (a random input that is not measured), and an ^ are, respectively, 

the output and input measurement noise. The measured response is y, and the 
measured input is "z. The noise sources w, t/q , and are assumed to be 

unknown random variables with zero mean. The transfer function H(w) is a 
matrix that defines the response to input at frequency to, that is, for 

u = u 0 e : ‘- a)t: the response is y = H(aj)uoe^ a)t . Then, 

H(co) = ^ = C(ia> - A) -1 B 


(For the applications in this report, only a single input and output are used, 
so the transfer function is a single function, not a matrix.) Let S^k) be 
the cross spectrum between z and y from the kth measured time-series of data, 
and s(k) an d s^k) the corresponding autospectra. Then, an estimator of the 
transfer function H(w) is given by the ratio of the averaged cross spectrum 
and input autospectrum 
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(ref. 12). 


H(a>) 


= 1EL = 

S 7. 


K 



k-1 


An analysis of the statistics of the estimator H is available in the 
literature (refs. 12 and 13). The variance of the spectra estimators S Z y and 
S z has been discussed previously for spectral analysis. The statistics of 
the ratio H are different, however, due to the correlation between the two 
signals. Bias errors in H are usually much less than the random errors (as 
long as the frequency resolution is chosen properly), so they will not be dis- 
cussed. The analysis of the variance of H gives the following result 
(ref. 12). Let a be the probability that the magnitude of the difference 
between the estimate H(oo) and the true value H(oj) is greater than the func- 
tion r(co). Thus, with confidence level (1 - a) , H lies with a circle of 
radius r about H, | H — H | < r, where 


-2 = 


K - 1 


F 2,2(K-l);a (1 


^zy) 


Sy 


Here, K is the number of averages; y Z y is the estimate of the input-output 
coherence function 



3 zy 1 


S z Sy 


and F 2 2 (K- 1 ) • a is t * ae va l ue °f the F-distribut ion with 2 and 2 (K - 1) 
degrees of freedom at the probability level 1 - a (i.e., the probability that 
F > F 2 2 (K_i); a is e 9 ua l to see ref. 10). For large K, the 

F-distribut ion may be approximated by F 2 2(K-l);a = _ ^ n a> Hence, the nor- 
malized error of the estimate of the transfer function (note [h| 2 = Sy/S z ) is 


s 


2 


-£n a 
K - 1 



) 


For a confidence level of 95 percent (a = 0.05), F = 3.1; at the 97.5 percent 
and 99 percent levels, F = 3.8 and 4.8, respectively. Assuming a Gaussian 
distribution, one standard deviation corresponds to a confidence level of 
68 percent, which gives F = 1.16. Hence, except for the effect of Y z y, the 
normalized error would be e 2 = 1/K, which is the same result as for the 
variance of the spectra alone. 

The coherence function y zy is a measure of how well y and z are 
correlated. With no noise in the system, Y Z y would be exactly equal to 1. 
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So, 1 - y__ 7 is a measure of the effect of noise on the calculation of the 
response due to a particular input. For the system described above, this 
factor is approximately 


1 - Y 


zy 


S z S y 


J zy I 


S z S y 


= s u ([h w 1 2 s w + s Vo ) + s vi s y 
(s u + s Vi )s y 




The first term is the ratio of the system response due to the process noise 
and the total response. The last two terms are the ratios of the measurement 
noise to the signal for the output and input, respectively. When y Z y is 
near unity, which is accomplished by having a high signal-to-noise ratio in 
the frequency range of interest, the normalized error of the transfer function 
estimate is reduced much more than with averaging alone. Thus, the number of 
averages, and so the total measurement time, needed to obtain a specified 
error level is reduced by correlating the input and output measurements. 


The transfer function is thus an accurate and efficient measure of system 
response. The experimenter controls the error level through the number of 
averages, and only the properties of the system itself are obtained, since 
both the input and output are measured. If the input and output are corre- 
lated (so the coherence function y Z y is near unity) , then the data required 
to achieve a specified level of accuracy is reduced. (In the experiment dis- 
cussed later in this report, the amount of data measured would have to be 
increased by a factor of 10 or 15 to achieve, by averaging alone, the accuracy 
obtained with the transfer function measurement.) Moreover, with the transfer 
function method, a quantitative estimate of the error level is available 
directly from the measured data 
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Hence, the dynamic stability measurement technique developed in this report is 
based on the system transfer function. 
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DAMPING RATIO FROM THE TRANSFER FUNCTION 


After measuring the transfer function, the experimenter must extract from 
H(oj) the parameters defining the system dynamics, in particular the frequency 
and damping ratio of the low-damped modes- A method is required to estimate 
accurately u) n and £ from the data in the vicinity of a resonant peak. This 
task is not easy because of the rapid variations in the transfer function mag- 
nitude and phase near the resonance, and because of the noisy transfer func- 
tion data. A number of least-squared-error curve-fit techniques was 
considered (fitting the data to a circle on the phase plane, or to a parabola 
on the Re H” 1 vs. Im H” 1 plane; or fitting the magnitude of the transfer 
function |h| vs. to) . Least-squared-error techniques generally do not work 
well with noisy data, and so the attempts to develop an algorithm along these 
lines were not successful. Sometimes a good damping coefficient estimate 
could be obtained, but usually the results for the- damping ratio were poor, 
even with clean data. 

The method adopted for analyzing the transfer function data is based on 
integrating H(u)) through a resonant peak. The use of integrals of the data 
reduces the sensitivity to noise in H. Consider the transfer function for a 
s ingle- degree-of- freedom, second-order system: 

= _1 = . _ JL 

m(o) n 2 - a) 2 ) + i C(jo m(o) n 2 - co 2 + i 2£od n u)) 

natural frequency, £ the damping ratio, C the damping 
m the modal mass. Now 


O) 2 do) _ 7T 

(oj 2 - o) n 2 ) 2 + (2£aj n u)) 2 ^Ca) n 

so it follows that the damping ratio is given by 



H 


where a) n is the 
coefficient, and 


C 


i U” im Ht ° dw ) 

n /" |H| 2 to 2 dw 


7TC0- 


and 

Im Ha) da) 

/ 0 " IhI 2 ^ 2 du 

The parameters o) n and £ define the eigenvalue of a mode, which is a property 
of the system as a whole, not of the particular input-output pair considered, 
while the modal mass and damping coefficient are dimensional parameters that 
scale with y/z. 
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In practice, the system has many degrees of freedom; in any case, the 
data are available only over a finite frequency range. Thus,, the integration 
is really performed over a finite band around a particular peak: 

03 = (1 - Af)oi n to a) = (1 + Af)u> n , where the bandwidth Af is a parameter of 
the algorithm. Then we have _ , 

TfO+Af^n 

_±_ U U-Af)ti»n 

0 7rw n |»(l+Af)a) n 

J (l-Af)o) n 

where the subscript 0 refers to the truncated integral. 

In the vicinity of a low-damped resonance, the principal effect of other 
modes is an overall phase shift. The expression above for the damping ratio 
was derived for a transfer function having a phase at the peak (co = 03p) of 


|h| 2 0) 2 du 


Im Ho) doj 




<H = -tan 


24 o)pO) n 

2 2 
w n - w p 


= -90° 


To apply this expression for £ , it is therefore necessary to rotate the 
measured transfer function on the phase plane until the phase at the resonant 
peak is -90° (see ref. 6). This is accomplished by introducing the following 
correction when integrating Im H 


H c = H 


-i Hp 
I Hr, I 


where Hp is the value of the transfer function at the peak (it is not 
necessary to correct the integral of |h| 2 , of course). 


Correction for Truncated Integration 


The effect of the finite limits of integration is that £ 0 is 1 GSS than 
the true damping ratio. To correct for the truncation of the integration, let 
£ = K£ 0 ; then 


K = 


da, 

J 0 D 


J 


( l+Af )a> n w 2 
( 1-Af )u) n D 


da) 


7T 

J ( l+Af )wn ^2 
( 1-Af )w n D 


da) 


where D = (a) n 2 - a) 2 ) 2 + (2£a) n a)) 2 . The integral in K can be evaluated 
analytically. Figure 1 shows the function K, obtained from calculations for 
various values of £ and Af; K is essentially a function of £/Af alone (or 
equivalently, Co/Af). A good approximation, which can be easily implemented, 
is 
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K = 


1 


2 io 

Af 



(see fig. 1). The first term is based on an evaluation of K for small £/Af 
the second term is added to account for the variation up to about £/Af = 2. 
The expression for the damping coefficient C does not require correction for 
the finite limits of integration. 


Numerical Integration 


The transfer function data are available at discrete frequencies, with 
separation Aw. Evaluating the integrals numerically gives then 


i 


(l+Af)oi) n 


( . . • ) du> 


(1-Af )u) n 


»£<••• 


) A03 


where the summation extends over the desired frequency range. The trapezoidal 
rule 


W n 


Aw 


Aw 


g do) = ^ (g 0 + §1 ) - g"(0 , <o 0 < 5 < u»i 


applied to Cq = (1/tt) 4£ 2 w n J (w 2 / D) dw, gives the integration error 


ACn » -M 


Aw 3 / 4C 2 w n ^ CJ 2 
12 


Now, the number of integration steps is M = 2Afw n /Aw; and the integrand 
curvature is (w 2 /D) m - -l/(2^ Lf w n i+ ) at the peak. Thus 




_ ( Af \ / Aw N 

The required frequency resolution for accurate integration is then 

Aw _ / 3tt 

Cw n ’ / Af ^0 


With Af = 0.1 and A£q “ 0.001 (a very good bound on the accuracy in calcu- 
lating £q), the requirement is Aw/£w n = 0.3. This is expected to be a con- 
servative estimate, since the curvature of the spectrum at the peak is used 
over the entire range to obtain ACq. Indeed, the results of numerical simu- 
lations of the algorithm indicate that Aw/£w n <.0.7 is adequate for negli- 
gible bias in C due to the integration error. Recall that the criterion for 
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negligible bias in measuring the spectra has exactly the same form, 

Aw/£w n constant. In general, the spectrum estimate bias is the stricter 
criterion, so satisfying it ensures the accuracy of the numerical integration 
in the damping ratio algorithm. 


Natural Frequency 

The other parameter required to define the eigenvalue is the natural 
frequency w n . The data available is the frequency w p at the r esonant peak. 
The natural frequency U) n is somewhat higher than w p (a) n = w p /A - 2£ 2 ) ; 
moreover, the true peak frequency is not what is available, since the transfer 
function is known only at points Aw apart. Based on the approximation 
|H| 2 = {2w n 2 mU 2 + (w/w n - l) 2 ]} 1 near the resonance, a corrected estimate 
of the natural frequency is 

Aw Hp 2 (Hr 2 - H l 2 ) 

Wn = t0p 2 V(V + H L F r- 2H R 2 H l 2 

where Hp = | H j at w = Wp, and Hj^p = J H | at w = Wp ± Aw. Numerical 
simulations of the damping ratio algorithm show this expression gives w n 
quite well. 


Alternative Expressions 

If only the magnitude of the transfer function is available, the damping 
ratio and damping coefficient may be obtained from 

= 1 SI |h[ 2 w 2 dw 

™n 3 |H P | 2 

and 


1 Hp | w n 

where | Hp | is the value at the resonance peak. The ideal transfer function 
also gives C = -Im H/[h[ 2 w for all frequencies (assuming negligible phase 
shift due to other modes). These expressions using the transfer function 
magnitude alone are also useful when the phase is distorted by other modes. 

Often, accelerometer measurements are used, rather than displacement or 
strain gage measurements. The damping ratio calculations may be applied to 
such data by simply using H a = -w 2 H, so 
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5 = 


as; 

T0) n r“ 

Jo 

and similarly for the other expressions. 


Im H a 


0) 

Ha| 2 


do) 


do) 


Damping Ratio Algorithm 


In summary, the damping ratio and natural frequency are obtained from the 
transfer function H(u)) as follows. The frequency u)p of a resonant peak is 
identified by examining the transfer function. Then, u) n and £ are calculated 
by the following expressions 


= 


Aw H P 2 (Hr2 - H L 2 ) 

“p + 2 ij_ 2 /'u_2 


Hp 2 (Hr 2 + H l 2 ) - 2 H r 2 H l 2 


■ h 


0.6366 

Af 


Cn + 


0.6817 


^ A -P / O ^ 


TCn'l 


Cn 


_ 1 _ [J(j-Af ^n Im HcU) dM ] 

I h ' 2 “ 2 


H c 


H * 


-i H p 

|H P | 


Similar expressions are given above to calculate the damping coefficient C, 
or to use the magnitude of the transfer function alone. Numerical simulations 
of the algorithm indicate that the calculations are not very sensitive to the 
choice of the integration bandwidth Af. Typically, a value Af = 0.1 is 
used, or less if necessary to avoid a neighboring resonance. (Note that the 
ratio of the integration bandwidth to the 1/2-power bandwidth equals Af /£.) 
Also, the correction for the finite limits of integration is an approximation 
valid up to about £/Af = 2, as shown in figure 1. The use of calculations 
based on the magnitude of the transfer function alone is recommended when a 
large distortion of the phase is due to other modes, to the extent that no 
circle on the phase plane exists, even approximately. This is often necessary 
for the higher modes of a system due to the increased modal density. Other- 
wise, the integral of the comp lex- valued transfer function should always be 
used, since it is less sensitive to noise in the data (particularly in | Hp | ) . 
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TRANSFER FUNCTION MEASUREMENT 


A typical test setup for the measurement of a system transfer function is 
shown in figure 2. A random signal generator provides the control input. The 
signal goes through a bandpass or lowpass filter so that the input energy is 
concentrated in the frequency range of interest. After establishing the test 
conditions and setting the signal generator parameters, the input is intro- 
duced into the system by increasing the potentiometer setting from zero, while 
monitoring the input and response. The input level finally established 
depends on what is required for a good signal-to-noise ratio, on the control 
actuator limits, and on the limitations of the system response. The potentiom- 
eter also quickly but smoothly removes the excitation, if necessary, such as 
for safety. Having established the excitation of the system, the online data 
analysis proceeds, beginning with the measurement of the transfer function 
between the input and some response. (The input signal may also be obtained 
from a transducer in the system, such as a servo position or load cell. In 
that case, the transfer function does not include the response of the input 
mechanism. The original input signal is used mainly when such an input meas- 
urement is not available, or when it shows a high level of measurement noise.) 
After the transfer function is measured, the excitation is removed by turning 
down the potentiometer. Then, the system parameters are calculated from the 
transfer function, perhaps while other testing takes place at this condition. 
After the stability level is obtained, and any significant trends assessed, 
the test condition may be changed. 

The transfer function may be measured in many ways. A digital computer 
offers great flexibility in obtaining the response data, and is essential for 
implementing the algorithm for extracting the system parameters from the data, 
even an algorithm as simple as the one used here. (A special purpose analog 
or hybrid computer could be constructed to do the job, but with considerable 
loss in flexibility.) A practical approach is to use a minicomputer system 
with specialized hardware devices, such as FFT and array processors. 

Figure 3 outlines the procedure for measuring the transfer function and 
calculating the modal frequency and damping according to the algorithm devel- 
oped in this report. While the random excitation of the system takes place, 
the selected input and output signals are sampled and digitized at a rate of 
r /sec. N samples are collected for each channel. The discrete Fourier 
transform is then applied to each of the two blocks of data. It is also con- 
venient at this point to convert the signals to engineering units. From the 
Fourier transforms of the input and output are obtained the autospectra and 
cross spectrum of the kth time series of data. The steps of data collection, 
Fourier transform, and spectra calculation are repeated K times, and the 
spectra averaged over the K records. Finally, the transfer function is 
obtained from the ratio of the averaged cross spectrum and the input auto- 
spectrum, and the error estimate e is calculated according to the expression 
given earlier. At this point, system excitation may be removed. 

Data analysis commences by examining the transfer function and designa- 
ting the frequency for a particular resonant peak. Then, the natural fre- 
quency and damping ratio are calculated, using the algorithm developed above 
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(the damping coefficient and modal mass may also be calculated if required). 
This analysis is repeated for all modes observable in the transfer function, 
and the resulting modal parameters are recorded. Finally, the transfer func- 
tion, autospectra, and normalized error may be examined and recorded as 
desired. 

The parameters defining the transfer function measurement are the sample 
rate r, the sample size N, and the number of averages K. The discrete 
spectra then have N/2 lines, with a maximum frequency (% ax = r/2 Hz (the 
Nyquist frequency). Thus, the frequency resolution is Acu = r/N Hz. The time 
required for a single record is Tq = N/r = l/Aw, and the total test time is 
T = K/A a). The considerations governing the choice of r, N, and K will now 
be discussed (see also refs. 12 and 13). 

The sample rate r should be about 2.5 times the bandwidth of interest. 
The last 20 percent of the resulting transfer function is therefore dis- 
regarded to allow for aliasing. Anti-aliasing filters, with a sharp cutoff at 
%ax ~ r /2, should also be used. Generally, the modes of interest should 
occur roughly in the middle of the spectrum. Data near the maximum frequency 
are suspect because of aliasing; data near zero frequency are suspect due to 
bias error. A system with many modes over a wide frequency range can be 
analyzed in successive runs with different bandwidths and frequency 
resolutions . 

The number of samples N equals r/Aa), so the choice of N is deter- 
mined by the required frequency resolution Au). The primary consideration in 
choosing Am is the bias error in the transfer function. The preceding dis- 
cussion showed that Agd = (1/2) £m n to £m n is required; the lower limit is 
usually conservative, so a value near the upper limit is typical. With a 
resonance in the middle of the spectrum (r = 4m n ) , the sample size should be 
N = r/Aoi = (4 to 8)/£, or N around 400 for a low-damped mode. Typically, 

N = 256 or 512 gives the required frequency resolution. In practice, the 
required sample size N should be established by trying two or three values 
at the beginning of the test, until a value for N is found that gives no 
noticeable bias error in H(rn) (the resonant peak amplitudes are the primary 
criterion) . 

The final parameter chosen is the number of records K, or equivalently 
the total data collection time T. In general, as many averages as possible 
are desired for the smallest error level; if necessary, the number of samples 
N could be reduced to shorten the test time. The averages required should be 
established by trying several values early in the test, and examining the 
error level as indicated by the normalized error e, noise in H(m), and vari- 
ations in the resulting parameters (especially the damping ratio). For a full- 
scale aeroelastic system, with low fundamental frequencies, the number of 
averages used is normally a compromise between what one would like for very 
low error and the requirements of an acceptable total test time. The more 
information that can be obtained from the data, the greater the justification 
for a longer test time. 
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EXAMPLE OF APPLICATION 


To demonstrate the application, the transfer function dynamic stability 
measurement technique was used in the flutter testing of a proprotor and 
cantilever wing model in a wind tunnel (fig- 4). The three-bladed, gimballed 
stif f-inplane rotor had a diameter of 0.86 m, and was operated at a speed of 
U = 22.67 Hz (1360 rpm) . The model was unpowered and hence was operating in a 
windmilling condition. The model was tested in the MIT Wright Brothers Wind 
Tunnel in axial flow over a speed range of V/UR = 0.34 to 0.69 (where V was 
the forward speed and UR was the rotor tip speed). The rotor controls con- 
sisted of collective and cyclic pitch actuators. The pylon also had a shaker 
vane (see fig. 4) that was used to excite the wing bending; the shaker vane 
was equivalent in effect to a wing flaperon. The fundamental wing frequencies 
were approximately 8.2 Hz (0.36/rev) for vertical bending, and 12.5 Hz 
(0.55/rev) for chordwise bending. ("Vertical" bending refers to the wing 
motion on the airplane, not to the orientation used in the wind-tunnel test.) 
Basically, the shaker vane was used to excite the wing vertical bending mode, 
and rotor collective pitch was used to excite the wing chordwise bending mode. 
Rotor cyclic excitation was also tried. Reference 14 gives complete details 
of the model. 

The proprotor and cantilever wing configuration is characterized by an 
aeroelastic instability at high inflow ratio due to the coupled motion of the 
rotor and wing. Thus, with this model, it was necessary to measure the fre- 
quency and damping of the fundamental wing modes, vertical and chordwise bend- 
ing, as a function of tunnel speed. What was desired in such a test was 
sufficient accuracy to detect trends in the stability so that the flutter 
boundary could be avoided, and enough information to allow correlation with 
stability predictions . Figure 5 shows typical results of the transient decay 
technique for this model. A scatter in the damping ratio data of A£ = 0.01 
or more was characteristic, so little useful information could be obtained 
with this method. The difficulty, which is typical of wind-tunnel testing, 
was that the response to tunnel turbulence obscured the transient, especially 
at high speed where the flutter measurement was most critical. Often there 
are more degrees of freedom involved than with this model, which makes it more 
difficult to produce and interpret the transient decay traces. Thus, the 
transfer function technique was used with this model to control the error 
level and thus obtain sufficiently accurate measurements of the system damping. 

Figure 6 shows the time history and autospectrum of the input signal used 
to measure the system transfer function. The input was white noise sent 
through a 3 to 12 Hz bandpass filter. Figure 7 gives the corresponding rms 
levels of the shaker vane and collective pitch used to excite the system. The 
level used was determined by the requirements for a significant response level, 
and the limitations on blade and wing Joads . The shader vane rms excitation 
was 1.5 to 2 deg, and the collective pitch was 0.25 to 0.75 deg; cyclic pitch 
was also used, at a level of about 0.25 deg (the cyclic pitch input was 
limited by an infinite-life blade load restriction). The one-half peak-to- 
peak level was around 3 to 4 times the rms value. These various control 
inputs were applied not simultaneously, but in sequential test points. Fig- 
ure 8 shows the resulting time histories of the wing and rotor motion. 
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comparing the cases with and without the external excitation (at V/ftR = 0.69). 
Figure 9 gives the corresponding rms levels of the wing vertical and chordwise 
bending response (or one-half peak-to-peak, which is 3 or 4 times the rms 
level; only the relative values with and without the input are of concern here). 
The wing response with external excitation, used in the measurement of the 
transfer function, was 2 or 3 times the response due to the existing tunnel 
turbulence (the process noise). Generally, the actual control position was 
used as the input measurement but, for a couple of the collective excitation 
points, it was necessary to use the original signal to the actuator. Fig- 
ure 10 shows the transfer functions of the rotor and shaker vane control actu- 
ators (the response magnitude has a linear scale) . The rotor control actua- 
tors had a response flat enough that it need not be considered in, obtaining 
the system transfer function. 

The model was tested at four speeds by successively exciting the shaker 
vane, rotor collective, and rotor cyclic pitch controls. The input (control 
position) and output (wing bending) signals were sent through 25 Hz low-pass, 
anti-aliasing filters. The signals were sampled at a rate of r = 51.2 Hz, 
with N = 512 samples collected for a single record. The spectra were aver- 
aged over K = 17 to 20 records. Thus, the spectra bandwidth was 
(%ax = 25.6 Hz; the frequency resolution, Aw - 0.1 Hz; the time for a single 
record, Tq = 10 sec; and the total data collection time, T = 170 to 200 sec. 

The damping ratio calculation used Af = 0.1. Figures 11 to 13 present the 
measured transfer functions at the four speeds, for, respectively, wing verti- 
cal bending response to shaker vane, chordwise bending response to collective, 
and vertical bending response to cyclic. The vertical scale of the transfer 
function magnitude is linear. The resonances of the fundamental wing modes 
are clearly evident. The scales of the transfer function magnitude plots are 
not the same in a given figure; the dimensional magnitude is not relevant for 
the present purposes anyway. Figure 14 gives examples of the transfer func- 
tion on the phase plane (Re H vs. Im H) ; the circle is characteristic of a low- 
damped resonance. Figure 15 shows the transfer function error estimates 
obtained from the data. The number of averages K = 17 to 20 gives a basic 
error level of e = 0.25. The lower values of e observed in figure 15 are 
the result of correlation between the input and output; the reduction in e 
occurs where there is a significant input level. To obtain this error level 
by averaging alone would require 10 or 15 times the data collected here, which 
would require a prohibitively long test time, even for a small-scale model. 

The efficiency of the transfer function measurement is clear. 

Finally, figure 16 presents the frequency and damping of the two wing 
modes obtained from the transfer function measurement by the algorithm devel- 
oped in this report. For the wing vertical bending, two measurements of the 
damping are available: from the shaker vane and rotor cyclic pitch excitation. 

At high speed, the shaker vane results are considered more accurate, as indi- 
cated by the line connecting the data in figure 16. This choice is based on 
the transfer function error levels found at high speed, as shown in figure 15. 
The cyclic transfer function shows a high error (due to limitations on the 
cyclic input level), while at low speed it shows an effect of the input-output 
correlation similar to that observed in the shaker vane transfer function. 
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This is an excellent example showing how the transfer function method cannot 
only control the error level but also directly assess it, and thus allow the 
experimenter to determine directly the most accurate measurements. 


CONCLUDING REMARKS 


The error statistics of several dynamic stability measurement techniques 
have been reviewed. The value of the transfer function measurement lies in 
the capability it gives the experimenter to control the error level, through 
averaging of the data and correlation between the input and output. Thus, a 
satisfactory error level may be achieved with a practical total test time. 

The characteristics of the system alone are obtained, independent of the input 
statistics (assuming an adequate bandwidth of the input spectrum). Moreover, 
the measured data give directly a quantitative estimate of the normalized 
error in the transfer function, which is extremely important in assessing the 
quality of the data. Thus, the transfer function is the basis for an accurate 
and efficient measurement of the system characteristics. An algorithm for 
obtaining the modal frequency and damping from the transfer function has been 
developed. The development of the test technique concluded by summarizing the 
rules for selecting the parameters involved in the transfer function 
measurement. 

The method developed here has two principal limitations. First, only a 
single- input , single-output transfer function was considered. Second, the 
data reduction algorithm is based on the response of a single-degree-of- 
freedom, second-order system; so, for multimode systems it is generally 
limited to analyzing the fairly low-damped modes. That is the most important 
information for flutter testing, but still the algorithm is only approximate, 
and the resulting information about the system is not complete. These limita- 
tions are fairly common in flutter-testing techniques. Therefore, a further 
development of the data reduction procedures would be desirable. The transfer 
function would still remain the basis of the stability measurement, providing 
the fundamental error control. The analysis procedure should handle multi- 
output, multi- input data; and it should completely and accurately identify all 
the poles and zeros of the system. Such a parameter identification algorithm 
would require more computational flexibility for the data analysis than the 
algorithm in this report. The major difficulties lie in developing an algo- 
rithm that accurately and efficiently obtains the system parameters from noisy 
data and that is still applicable for on-line analysis in flutter testing. 

The assumptions inherent in the use of a transfer function representation 
are that the system is linear and t ime- invar iant . These assumptions are quite 
satisfactory in practice, but time-invariance leads to further problems with 
rotorcraft testing. To have a time-invariant (constant coefficient) system 
often requires nonrotating degrees of freedom representing the rotor motion, 
which in turn requires that the motion of all the rotor blades be measured. 

Such complete measurements are seldom available; usually, only one blade is 
instrumented. What is required in such cases is a filter (i.e., an output 
error estimation method) that estimates the nonrotating states of the rotor 
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from measurements of the individual blades in the rotating frame, especially 
with data from only one or two of the blades. Then, the transfer function 
methods may be applied. 
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APPENDIX 


ERROR ANALYSIS OF THE METHOD OF RANDOM 
DECREMENT SIGNATURES 


The method of random decrement signatures (randomdec) is a procedure 
developed by Cole (ref. 15) for analyzing the response of a linear system to 
random disturbances. The procedure is designed to estimate the transient 
response of the system. This appendix analyzes the error statistics of the 
random decrement signatures. The randomdec estimator y is defined as the 
ensemble average of the response to existing random disturbances 

y (t) = rE ^ (t) 

k=l 

As usual, the ergodic hypothesis is used, and the random decrement signature 
is obtained by averaging records that are sequential in time 

5(t) ‘IEK t + t « k> ) 

k=l 

The records are assumed to be uncorrelated, which is generally true if there 
is no overlap (ref. 15). The key to the procedure is to select each ensemble 
with the same initial conditions by an appropriate triggering method. 

Then, all records have an identical transient due to the initial conditions, 
while the subsequent noise averages out. 

Consider again the linear system x = Ax + Bu, and response measurement 
y = Cx + v, excited by the unknown random disturbance $, with measurement 
noise v. Assume that u and v have zero mean. The solution using the ini- 
tial conditions at time to is 

— >- n A ( t t o )^ / , \ , n f t A ( t t ) j . ->■ 
y = Ce u x(t 0 ) + CJ e Bu dx + v 

to 

The observation y(t) is thus a random process, depending on the random vari- 
ables u, v, and x(t 0 ). The expected value of the random decrement signature 
is then 

£y = = Ce ACt_t % 0 

k=l 

Here, the initial condition ^(t 0 ) is a stationary random variable with 
expected value xg = Ex(tg) . Hence, y is an unbiased estimator of the system 
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transient response. It also follows that y depends on the trigger procedure, 
which determines x 0 . Clearly, the frequency content of the random disturb- 
ance u must be sufficient to excite the modes of interest, for xq must be 
selected from the steady-state response to u. Moreover, the trigger process 
must be such that there are nonzero initial conditions for the modes of inter- 
est. The usual procedure is to trigger the start of the ensemble on a fixed, 
nonzero level of the response y, with either positive or negative slope. 

Cole (ref. 15) suggests a trigger level of y = 1.0 to 1.6 Oy (where Oy is 
the rms level of the response) . 

The variance of the random decrement signature y is 


Vy = Ey 2 - (Ey) 2 

= C e A(t - t 0 ) V X() e AT(t " t 0 ) C T 

+ ^ jcj^ e A ^ t:_T ^BQB T e AT< - t_T - > dxC T + R j 


Here, it has been assumed that u and v are uncorrelated; that the measure- 
ment noise covariance is R; and that the excitation u is white noise with 
Eu(ti)u t (t 2 ) ~ Qd(xi - t 2 ) . The ensembles are uncorrelated, except for the 
initial conditions. The triggering process determines x(to), with expected 
value xq - E x ( t o ) , and variance Vxq = Ei(tg)x T (to) - XqxJ. (The steady- 
state response of the system to the disturbance u has zero mean and variance 
a x 2 = Ex(t)x^(t), which is not the same as xq and Vxq due to the procedure 
for selecting the initial time tg . ) The second term in Vy is the same 
error as found for the transient decay, divided by the number of averages K. 
It follows that for a signature with a transient that decays to a fraction f 
of its initial value, this term gives a maximum error approximately 


e 


2 



(cf. the analysis above for transient decay). So, the ensemble averaging 
reduces the error in the transient due to the response of the system to the 
random disturbance u subsequent to the time tg where the initial condi- 
tions are determined. As usual, the error is inversely proportional to the 
square root of the number of averages. This error source can be made as small 
as required by increasing the number of averages. (Here, the experimenter 
does not have control of a x /xg, since xo must be obtained from the time 
series with rms level a x ; generally, a x /xg = 0.5 to 1.) 

The remaining term in Vy depends upon the variance of the initial con- 
ditions Vxq . Both this term and (Ey) 2 decrease with time at the same rate, 
so the normalized error in y always has about the same value, equal to the 
normalized error in x(tg) 
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2 Vy _ Vx 0 
E (Ey) 2 x Q 2 

The ensemble averaging does not affect this term, so an accurate response 
measurement depends upon keeping the variance of x(to) as low as possible. 

The triggering process is required to select to based on the measurement of 
y, such that x = xq . An error in y results unless exactly the same initial 
conditions are produced by every triggering event. The sources of the vari- 
ance Vx 0 include: variation in the trigger mechanism or algorithm; several 

modes observable in y; and measurement noise in y. The measurement noise is 
probably the most important problem; with a trigger based on y, the measure- 
ment noise covariance must give a lower bound on Vxq . With high measurement 
noise, the experimenter can do little to control the error in the random decre- 
ment signature estimate of the transient response. Increasing the number of 
averages does not influence this error; the trigger level cannot be increased 
much without increasing the correlation of the records (ref. 15); and it is 
not often possible to increase the excitation u relative to the measurement 
noise 

Finally, consider the total amount of data required. The time for a 
single record, with decay to a fraction f of the initial value, is 
To = -ftn f / (2 tt C to n ) - Then, assuming that the Vxq error is negligible, the 
total sample time required is 


T = KT 0 = 


constant 

Ccu n e 2 


This is the same result as for spectral analysis (allowing for a somewhat dif- 
ferent meaning of e; see also ref. 15). It is the general result for the 
amount of data required to estimate the statistics of a random signal. As for 
spectral analysis, for a given T it is necessary to compromise between fre- 
quency resolution and the number of averages. Typical applications of the 
random decrement signature procedure require K = 2000 to 8000 averages for 
satisfactory accuracy (refs. 15 and 16). Data processing for the randomdec 
method follows the lines discussed for obtaining the damping from a transient 
decay trace, with improved accuracy due to the ensemble averaging. Analysis 
in the time domain is rather more sensitive to noisy data than analysis in the 
frequency domain, but it has been successfully implemented in practice. 

As with autospectral analysis, there are two basic difficulties with the 
method of random decrement signatures. First, there is some uncertainty 
whether only the properties of the aeroelastic system are being estimated 
(because only the system output is measured). In particular, the initial con- 
dition variance is an error source that can neither be controlled by averaging 
nor be estimated directly from the experimental data. Second, the total meas- 
urement time required to obtain accurate estimates of the system properties is 
likely to be very long. This large amount of data required to extract informa- 
tion accurately from the signal is the basic characteristic of the class of 
methods utilizing response measurements alone. For further discussion for the 
random decrement statistics, see reference 11. 
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Figure 3.- Outline of dynamic stability measurement procedure based on the 

transfer function. 
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Figure 5.- Examples of transient decay traces for the proprotor and wing model 
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7.- BMS levels of shaker vane and rotor collective pitch inputs used to 

excite the model. 
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Figure 8.- Time histories of wing and rotor response with no input (wind- 
tunnel turbulence only) and with external input to excite the model at 
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Figure 14.- Typical transfer functions on the phase plane for wing vertical 
bending response to shaker vane (q x /6 F ) and for wing chordwise bending 
response to collective pitch (q 2 /0 q ) at V/ftR = 0.69. 
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Figure 16.- Natural frequency and damping ratio of the wing vertical and chord- 
wise bending modes obtained from the transfer function measurements. 
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